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Abstract 

In this paper we compute the four-loop corrections to the QED photon self- 
energy H(Q 2 ) in the two limits of q = and Q 2 — > oo. These results are 
used to explicitly construct the conversion relations between the QED charge 
renormalized in on-shell(OS) and MS scheme. Using these relations and results 
of [T] we construct the momentum dependent part of the H{Q 2 , m, a) at large 
Q 2 at five loops in both MS and OS schemes. As a direct consequence we 
arrive at the full result for the QED /3-function in the OS scheme at five loops. 
These results are applied, in turn, to analytically evaluate a class of asymptotic 
contributions to the muon anomaly at five and six loops. 



1 Permanent address: Institute for Nuclear Research, Russian Academy of Sciences, Moscow 117312, 
Russia. 



1 Introduction 



Quantum electrodynamics (QED) is one of the best tested and established quantum field 
theories. The study of its input parameters and properties is thus a challenge for each 
theoretician and experimentalist. For example the anomalous magnetic moment of the 
muon a M has been measured with impressing accuracy at the level of 0.5 parts per mil- 
lion |2,3J: a™ p = 116592089(63) • 10~ n . Also from theory side the anomalous magnetic 
moment has been studied in great detail through the computation of higher order correc- 
tions. In general higher order corrections to a^ eo are classified into three classes: pure 
QED, electroweak and hadronic contributions, where sample diagrams are depicted in 
Fig.m 




Figure 1: Sample diagrams for the classification into: pure QED contributions (a), elec- 
troweak corrections (6i), (62), and hadronic contributions (c). 

Within this work we consider higher order corrections to the pure QED part. Recently 
the complete tenth order QED contribution has been determined numerically in Ref. [I]. 
We refer to the reviews [5H7] and references therein for a discussion of the lower order 
QED corrections, the electroweak and hadronic contributions. 

In QED diagrams with internal fermion loops arise starting from two-loop order, where 
the fermion type of the internal loop can be in general different from the external muon. In 
the case that the internal fermion loop consists of an electron, logarithmic contributions 
of the type ln(M M /M e ) arise, where is the mass of the muon and M e the mass of 
the electron. In view of the large mass ratio M ti jM e ~ 200 one can expect, that these 
logarithms play a dominant role. These logarithmically enhanced contributions arise on 
the one hand from the insertion of the electron vacuum polarization (eVP) into the first 
order muon vertex diagram, shown in Fig. [T](a), but they can on the other hand also 
appear through light-by- light (LBL) scattering diagrams. Examples for both diagram 
types are shown in Fig. [2j 




Figure 2: Example diagrams leading to dominant logarithmic contributions from electron 
vacuum polarization insertions and light-by- light scattering diagrams. 
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It has been demonstrated in Ref . [8] , that the asymptotic part of the anomalous magnetic 
moment of the muon a a ^ ymp , which contains these logarithmic contributions originating 
from the electron vacuum polarization function insertions, as well as the mass independent 
term can be obtained with the help of the electron vacuum polarization function in the 
asymptotic limit M e — > 0. This technique has been applied in order to derive a°^ ymp up to 
three-loop and four- loop order in Refs. [8HTT]. A similar technique has also been applied 
in Ref. [12]. 

At five-loop level all diagrams contributing to a l ^ eo can be decomposed into six super- 
sets, which can be further subdivided into 32 gauge invariant subsets [13] • These supersets 
have been computed numerically in Ref. [4]. For a few subsets also analytical results could 
be obtained in Refs. [T4j[T5]. 

We will compute analytically the asymptotic limit of the photon propagator with 
massive electron loops and photon exchanges at four-loop order in two renormalization 
schemes: the MS-one and the classical on-shell scheme. The calculation is then employed 
for the derivation of the MS-on-shell relation of the fine structure constant a at four- 
loop order. This relation allows a conversion of observables, conveniently computed in 
the MS scheme into the on-shell(OS) scheme or vice versa. These conversion relations 
combined with the recently published QED /3-function in the MS scheme at five loops [1] 
are then used to derive the complete five-loop contribution to the QED /3-function in 
the OS scheme as well as the momentum- dependent part of the polarization function in 
both the MS and OS scheme at five loops. These new four-loop and five-loop results are 
subsequently used to determine analytically a a ^ ymp at five loops as well as some genuinely 
six- loop contributions to a a ^ ymp . This will constitute a new analytical result for several 
gauge invariant subclasses of the five-loop QED muon anomaly and will serve as a check 
for some already known numerical results. 

The outline of the paper is as follows: in Section [2J we introduce our notations and 
conventions. In Section [3] we discuss the methods of calculation and present the results 
for the vacuum polarization function for small and large Q 2 in MS and on-shell scheme. 
In the next Section we use these results to derive the conversion formulas between the 
QED charge, renormalized in MS and on-shell schemes. In Section [5] the polarization 
operator at five loops is discussed. In the following Section [6] we derive the complete 
five-loop contribution to the QED /3-function in the OS scheme as well as the momentum- 
dependent part of the polarization function in both MS and OS schemes at five loops. 
In Sections [7] and [8] we analytically compute all asymptotic contributions to the muon 
anomaly at five loops and logarithmically enhanced ones at six loops respectively. Finally 
in Section M we close with a brief summary and our conclusions. 

2 Notation and Generalities 

We consider QED with N identical fermions with mass m. The MS renormalized photon 
propagator D R is related to the photon polarization function H(—q 2 //j, 2 ,m/fx,a) in the 
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standard way: 
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r \ 
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(/i) and /i are the running coupling constant, the fermion 



where a 

mass and normalization scale in MS scheme respectively. The symbol q denotes the 
external Minkowskian momentum of the polarization function. At four-loop order the 
appearing self-energy diagrams can be classified into four classes: singlet and non-singlet 
type diagrams, where the non-singlet ones can be again decomposed with respect to the 
number of inserted closed fermion loops. Some example diagrams are shown in Fig. |3j 




Figure 3: The four diagram classes arising in the computation of the four-loop vacuum 
polarization function are shown. Only one representative for each class is given. Class (s) 
are singlet type diagrams, whereas class (a), (b) and (c) are non-singlet type diagrams, 
which are decomposed into the number of inserted fermion loops. 



As is well-known the MS renormalized photon self-energy does not have any power-not 
suppressed terms of the type ln(m 2 ) in the limit of large Q 2 = — q 2 . They appear only at 
order C(m 4 /Q 4 ) of the large Q asymptotic expansion (see, e.g. Ref. [16J). This allows to 
define the asymptotic part of the polarization function II in the limit m — )■ as follows: 
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The running of the coupling constant is governed by the corresponding evolution 
equation of the form: 



fi 2 -^\na = (3(a). 



(3) 



The evolution equation for the polarization function can be directly obtained from the 
fundamental concept of the scheme-invariant charge [TTIIT8] . defined through 
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For the case of the asymptotic polarization function, the evolution equation assumes 
especially simple form: 
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Since recently the QED /3-function (3(a) is known to five-loop order pQ. For future refer- 
ence we provide the corresponding result belowQ 
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Traditionally in calculations of leptonic(£) anomalies ae the classical OS scheme is 
employed. In this subtraction scheme the QED is parametrized with two variables, namely, 
the fine structure constant a = a os , equal to the invariant charge (jl]) at zero momentum 
transfer, and the pole mass M of the fermion. The corresponding normalization condition 
for the OS polarization function q 2 /M 2 ,a) = Il os (— q 2 /M 2 , a) is: 



n(-g 2 = 0,M,a) = 0. 



(7) 



The asymptotic polarization function U asymp (Q 2 /M 2 , a) is obtained from I1(Q 2 /M 2 , a) by 
discarding in every order of perturbation theory all power suppressed terms in the limit 
of M — y 0. Finally, according to Ref. [8] one obtains the asymptotic part of the muon 
anomaly (related to the vacuum polarization diagrams) through the relation: 
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where 



4 Symp (Q 2 /M 2 ,a) 



i + n as y m p(Q 2 /M 2 ,a)' 
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3 The polarization function at four loops 

The asymptotic limit m — > requires the computation of massless propagators, whereas 
the limit g 2 — >■ leads to the evaluation of massive tadpole diagrams. The perturbative 
expansion of the vacuum polarization function in the fine structure constant in the two 
cases is conveniently defined by: 

n(g 2 ,m 2 = o,a) = En (i) (W (^-J (10) 

2 Note, please, that the function /3(a) is related to f3^ ED (A) as given in Eq. (4.7) of p] as follows: 

m = (pq ed (a)/a) u =s/(4w) . 
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n(Q 2 
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where £ m q = ln(/i /Q ) and l^m = ln(/^ /m ). In both cases the computation has been 
performed with FORM [I~9Tl2~l~] based programs. The four-loop QED diagrams contributing 
to the electron polarization function have been generated with the program QGRAF [22J. 

The four-loop massless propagators can be reduced to 28 master integrals. The re- 
duction of four-loop massless propagators has been done by evaluating sufficiently many 
terms of the 1/D expansion |23j of the corresponding coefficient functions [24J. The master 
integrals are known analytically from [25||26]. 

In the low Q 2 limit all appearing tadpole diagrams have been reduced to master 
integrals with the help of Laporta's algorithm [271128] . The arising polynomials in the 
space-time dimension d = 4 — 2 e have been simplified with the program FERMAT [29] . The 
remaining master integrals are known analytically to sufficient high order in e and have 
been taken from Refs. [30H38] . 

For the large Q 2 limit our result, renormalized completely in the MS scheme, reads: 
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with being the Riemann zeta-function defined by: 
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and £ m q = ln(^a). The symbol si is here and in the following equal one and serves only 
as a separator in order to display the singlet contributions, originating from diagrams of 
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the type as shown in Fig. |3]^s), separately. For large Q 2 the MS renormalized photon 
self-energy does not have any logarithms of the type ln(m 2 ). They appear only at order 
Oim^/Q 4 ). However, in the low Q 2 = case the renormalized vacuum polarization 

2 

function does contain logarithms of the type £^ m = ln(=r). In this limit our analytical 
result, renormalized in the MS scheme, reads: 
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where the polylogarithm function Li n (l/2) is defined by: 
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In the next step we convert the results of Eqs. ( I12T) - (TT5T) from the MS into the 
on-shell scheme in order to obtain the asymptotic limit of the polarization function 
n asymp (Q 2 , M 2 , a) in OS scheme. The scheme independence of the invariant charge di- 
rectly leads to the relation: 
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where it is understood that the running a on the r.h.s. is converted into the on-shell 
coupling constant a. The corresponding conversion relation is given in the next section HI 
In analogy to Eqs. (fTOl) and ( fill) the expansion in the on-shell fine structure constant 
a of n asymp (Q 2 , M 2 , a) is defined by: 
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The asymptotic limit of the polarization function in the OS scheme is given by {Imq 
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This result in OS scheme does not contain any logarithms of the renormalization scale [i 
anymore. 
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4 The QED charge renormalized in MS and on-shell 
scheme 

Let us define the conversion factor C a a, which converts the fine structure constant a in MS 
scheme into a in OS scheme: a = C a ^a. The conversion factor C Q „ has a perturbative 
expansion in a defined by: 



The expansion coefficients can be obtained by evaluating Eq. ( 122]) at Q 2 = and 
using Eq. ([7]). On the r.h.s. of this equation we insert the results of Eqs. f lTTj) . f lT7|) - 
( l20l) and expand in a up to four-loop order. In addition, one should use the relation to 
convert the MS-mass m(/z) to the on-shell mass M, which is known from Refs. |39j[l0] to 
three-loop order. For the case of QED it is given in Appendix [Al 
The different orders read: 
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with Cm = ln(/i 2 /M 2 ). On the other hand the inverse conversion factor Caa, which allows 
a conversion from the MS to the on-shell scheme a = Caa ot is useful as well. In analogy 
to Eq. (128p the perturbative expansion in a of Caa is defined as: 
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Proceeding like in the previous case and inverting the series obtained with the help of 
Eq. (J22J , we find: 
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5 The polarization function at five loops 

In this section we will use the recent calculation [1] of (3(a) — the QED /3-function in the 
MS scheme — at five loop order to find the Q-dependent part of the five-loop contribution 
to the asymptotic polarization function in both MS and OS schemes. 



5.1 MS scheme 

We start by transforming evolution equation (JSJ) for the MS-renormalized asymptotic 
photon polarization into the form: 
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Now by direct integration of the r.h.s. of Eq. ( 138)) one can easily construct the Q-dependent 
part of n asymp (%-,'aj in order a 5 : 
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5.2 OS scheme 

By inspecting Eq. f l22|) relating the photon polarization function normalized in the MS 
and OS schemes one observes that the Q-dependent part of the 0(a 5 ) OS-normalized 
asymptotic polarization function is recoverable from Eq. ([6]), four-loop conversion relations 
and five-loop asymptotic MS polarization function. The result is: 
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6 The QED /3-function in the OS scheme at five loops 



The OS /3-function, 
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describes the evolution of the asymptotic photon polarization function with respect to the 
change of the mass M, namely (see, e.g. [4Tj): 
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At the four-loop level /3(a) is known since long One could easily check that Eq. (HZ 
is indeed met by our four-loop result given in Eqs. f )23|) and (j2 
For our purposes it is useful to rewrite Eq. ( 142]) in the form 
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M dAP + ^ a)a d^ a 



/3(a) 



a 



(43) 



Eq. (143 p demonstrates that one could construct the next (L + 1) order contribution to (3 
provided one has the polarization function n asymp in L-loops (that is up to and including 
the terms of order a L ) and its all logarithmic (that is proportional to C-mq) terms of order 
a L+1 . Thus, we have at our disposal all ingredients to find /3 at five loops. The result 
reads: 
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(48) 



7 Contributions to the muon anomaly at five loops 

Having at hand the result for the asymptotic polarization function at four-loop order 
in QED, one can use it in combination with Eq. fl9]) and Eq. (JS} to determine the full 
asymptotic contribution (that is coming from electron vacuum polarization insertions in 
the limit M e /M M — > 0) to the muon anomaly at five- loops. The first of the six supersets 
of diagrams, which contribute at five-loop order to the muon anomaly is shown in Fig. HI 
It can be subdivided into ten gauge invariant subsets. 




1(a) 




Figure 4: The ten gauge invariant subsets contributing to the muon anomaly which 
originate from inserting the vacuum polarization up to four-loop order into the first order 
QED vertex of Fig. Ufa). For each diagram class only one typical representative is shown. 
Wavy lines denote photons(7), solid lines denote electrons(e) or muons (/i). The last 
five diagrams {1(f), 1(g), 1(h), I(i), I(j)} are non-factorizable insertions of the vacuum 
polarization function; the first five diagrams {1(a), 1(b), 1(c), 1(d), 1(e)} are factorizable 
ones. 



We define the perturbative expansion of the anomalous magnetic moment of the muon 
by: 



asymp 



E 



r asymp,(i) ( ^_ 



1,2,3 



The decomposition of a°f ymp into the ten subsets shown in Fig. H]is given by: 



a 



asymp, (5) 



FigM 



+ a (f)>{/} + a (5).{*+fc> + a (5),« + a (5),m 



(49) 



(50) 



The letter in the upper case curly brackets denotes from which contribution of Fig. |4] 
the considered term is coming. Inserting the results of Eqs. (l23l)- (|271) into Eq. (jBJ) and 
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performing the integration the results for the different contributions of Eq. (|50|) read: 
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(51) 

(52) 
(53) 

*) 

(54) 
(55) 

(56) 



(57) 



(58) 



(59) 



with 



ln(M M /M e ) and power suppressed terms of the order 0{M Z IM^) are neglected 
in the asymptotic limit. The coefficients for the logarithmic terms have been computed 
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analytically in Ref. [H], except for the complete ^ e -term in Eq. ( 155]) . For the factorizable 
insertions of the vacuum polarization function also the mass independent term has been 
determined analytically in Refs. [T41IT5] . except for the case /(e). All the analytical 
results of Refs. (THEE] for the Eqs. ( )56l) -(l5Tl) are in agreement with ours. Numerical 
results for these ten gauge invariant subsets of diagrams which are shown in Fig. H] have 
been reported in Refs. [4*2TI4"5] . A comparison between the asymptotic analytical formulas 
of Eqs. fl5T]) -f l56|) and the numerical results is shown in Tab. [TJ 



Subset 


analytical 


numerical 


Ref. 


num.-ana. 


1(a) 


20.1832 + 0(M e /M tl ) 


20.14293(23) 


[42j 




-0.04 


1(b) 


27.7188 + 0(M e /Mj) 


27.69038(30) 


m 




-0.03 


1(c) 


4.81759 + 0(M e /Mj) 


4.74212(14) 


m 




-0.08 


1(d) 


7.44918 + O(MJmJ) 


7.45173(101) 


m 




0.003 


1(e) 


-1.33141 + 0(M e /Mj) 


-1.20841(70) 


[42j 




0.12 


Hf) 


2.89019 + 0(M e /M il ) 


2.88598(9) 


[A2\ 




-0.004 


1(g) + 1(h) 


1.50112 + 0(M e /Mj) 


1.56070(64) 


m 




0.06 


m 


0.25237 + 0(M e /M l }j 


0.0871(59) 






-0.17 


m 


-1.21429 + 0(M e /M l j) 


-1.24726(12) 


[45j 




-0.03 



Table 1: The first column shows the different gauge invariant subsets of diagrams as 
defined in Fig. HI The second column contains the corresponding results of Eqs. (!5T!) - (l59|) 
evaluated numerically, where we have used for the mass ratio M^/M e = 206.7682843(52) 
|46j . This result is correct only up to power corrections in the small mass ratio M e /M^. 
The third column contains the numerical result obtained in Refs. [4"2~H4"5] . (Note that in 
Ref. [42J also a more precise value 7.45270(88) is given for 1(d) which was obtained using 
the exact sixth order spectral function.) The last column shows the difference between 
the numerical and asymptotic analytical results. The subsets {1(a), 1(b), 1(c), 1(d), 
1(e)} originate from Feynman diagrams with factorizable vacuum polarization insertions, 
whereas the subsets {1(f), 1(g), 1(h), I(i), I(j)} are non-factorizable (see Fig. H]). 



In general good agreement between asymptotic analytical and numerical results is found, 
except for the subsets 1(e) and I(i) where we only have poor agreement. The remaining 
differences should arise from corrections of the order 0(M e /M fl )^ Summing up all ten 
subsets one obtains a a ^ ymp ^ F - g] = 62.26675 + 0(M e /M^) which is, despite the small 

number of diagrams, sizeable (~ 8%) compared to the complete numerical result a® = 
753.29(1.04) of Ref. jl]. The reason for this is the logarithmic enhancement. Indeed, if we 
set to zero all Q-dependent terms in the asymptotic OS polarization function the result 
for a ™ ymp ^ r-n would be reduced to 2.52261. 

^ Fig. |4J 



! For a discussion of the differences for subset I(i) see also Section VIII of Ref. [44] 
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8 Contributions to the muon anomaly at six loops 

The five-loop contribution to Eq. (jHJ) 

dT mp (Q 2 /M" } a) = Y,d l (i MQ ) (- 



i>0 

can be schematically represented as follows: 

d 5 = 4 actr - n asymp - (5) + C 5 , 

where 



(60) 



d 



factr 



_ ^Jjasymp,(l)^ _|_ ^ ^-Qasymp,(l)~j j-jasymp,(2) _ g jjasymp,(l) ^-Qasymp,(2) x 

3 fn asymp, ( 1 ^ 2 n asymp ^ 3 - ) + 2 n asymp,( - 2 ^ n asymp,( - 3 ^ + 2 n asymp ^ 1 - ) n asymp ' ( ' 4 * ) (61) 



and C5 is a still unknown constant standing for Q-independent contributions missing in 
n asym P ,(5) ag Ascribed by Eq. (jlO]). 

The corresponding result for the asymptotic contribution to the muon anomaly is 
rather bulky so we provide it below in numerical form only: 



a asymp,(6) = 257.245 + — = 246.381 factr + 10.8647 + — . 



(62) 



where after the second equality sign we decompose the result into two pieces: the factoriz- 
able and the genuine six-loop terms. We observe that the six-loop asymptotic contribution 
to the muon anomaly is almost completely saturated by factorizable terms. One could 
hardly expect that the constant C5 might have any numerical relevance. Indeed, an analog 
of Eq. f[6"2"j) at the five-loop level looks as: 



ar ,m R (5) = 60^524 + Ci - 3.0287 = 5M374 ,^ + lm5 + 9 - 3.0287 



Thus, we believe that Eq. fl62|) with C 5 = presents quite a good prediction for a°f ymp ^ 6 \ 



Subset 


AT 4 


N 3 


N 2 


N 


siiV 3 


siiV 2 


Value 


7.15995 


7.15320 


-2.73813 


3.37488 


-7.22943 


3.14426 



Table 2: The numerical values are correct up to power corrections of the order O (M e /M^) 
and the contribution coming from the unknown constant C5 of Eq. ( 160]) . The symbol si 
labels singlet contributions, whereas the power of N denotes the number of closed fermion 
loops which arise in the corresponding diagrams. 

To further study the size of the factorizable and non-factorizable six-loop contribu- 
tions to the muon anomaly we subdivide them into several gauge invariant subsets in 
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complete analogy to the five-loop case shown in Fig. HI Let us start with the genuine, 
non-factorizable contributions which correspond to the term — n asymp ' 5 + C5 of Eq. (1601) 
and which we again subdivide according to the number of fermion loops arising in the cor- 
responding diagrams. The individual contributions to a " symp ' (6) = 10.8647 + CJ2 

r non— factr 

of Eq. ( 1621) are shown in numerical evaluated form in Tab. [2] and an example diagram for 
each of the six subsets of Tab. [2] is given in Fig. [51 




Figure 5: For each subset of Tab. |2]one example diagram is shown. It would be possible to 
subdivide these diagram classes into further subsets, however, we refrain from doing this 
for simplicity. The symbol si denotes a singlet diagram, whereas the power of N gives the 
number of closed fermion loops which arise in the corresponding diagrams. Wavy lines 
represent again photons whereas solid lines stand for electrons or muons. 



The factorizable contribution of Eq. 



a asym P ,(e) _ 246.381, which originates 

factr & 

from the terms in of Eq. (I6TT) can also be classified into several subsets: 



a 



asymp,(6) 



factr A* M M 

a (6),M +a (6),{M + a (6),{^} 

£t fj, £1 



+ 



a W.W + a (6),{e> + a (6),{/} 

fj, fj, fj, 

a (6) 1 W + a (6), Wj 



(64) 



where we have shown in Fig. [6] one typical contributing diagram for each term on the r.h.s. 
of Eq. (16~4"1) . The letter in the upper case curly brackets denotes again from which contri- 
bution of Fig. [6] the considered term is coming. The numerical values of the factorizable 
contributions are shown in Tab. [31 
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I 6 (a) l e (b) T 6 (c) I 6 (d) 




left leti) IbW ' leW 



Figure 6: The different diagram classes of the factorizable six-loop contribution. For each 
diagram class I§(o)-Iq(1) one typical representative is shown. 



Subset 


J 6 (a) 


h{b) 


h{c) 


h{d) 


h{e) 


W) 


Value 


58.1861 


101.150 


35.8953 


30.6997 


-5.12107 


9.68427 


Subset 


h{g) 


h{h) 




h{k) 


IS) 


Value 


-1.81937 


15.6754 


7.73424 


0.68395 


-6.38790 



Table 3: Numerical evaluation for the individual factorizable contributions of Eq. (|64|) . 
The numerical values are correct up to power corrections of the order O (M e /M^) and 
arise from the diagram classes shown in Fig. [61 

9 Summary and Conclusion 

We have computed analytically the vacuum polarization function at four-loop order in 
perturbative QED in the asymptotic small and large momentum limit. From the low 
energy limit we have derived the conversion factor of the fine structure constant between 
on-shell and MS scheme. These results have been used to derive an analytical expressions 
for the dominant and gauge invariant contributions to the muon anomaly originating 
from vacuum polarization function insertions at five-loop order. Numerical results for 
these contributions are already available in the literature and are up to power corrections 
mainly in agreement with our analytical ones. 

Using the result of the recently computed five-loop QED /3-function in the MS 
scheme of Ref. [1] we also determine the asymptotic momentum dependent part of the 
polarization function at five loops which in turn allows to calculate vacuum polarization 
type, asymptotic, leading contributions to the anomalous magnetic moment of the muon 
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at six loops. The five-loop momentum dependent part of the polarization function is also 
used to determine the five-loop QED /3-function in the on-shell scheme. 

All the calculations described here were done during 2007-2008. A short version of 
this work, dealing only with the N=l QED was reported at the 9th DESY Workshop on 
Elementary Particle Theory in the spring of 2008 [17] . 

Our results will be made available in computer readable form under the URL: 
http : //www-ttp . physik . uni-karlsruhe . de/Progdata/ttpl2/ttpl2-021 . 
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A MS— on-shell relation for fermion masses in QED 

The relation between the MS-mass m(/i) and the on-shell mass M for N identical fermions 
in QED is given by 1351130): 




1^ 9 1 1 

+ ^ N£, M + - i\ M + - N q M + - tt 2 ln(2) 



+ 



( 
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+ -iV7r 2 ln(2) 



vr 2 i m ln(2) - l - N vr 2 l» M ln(2) + 1 vr 2 ln 2 (2) 



^iV7r 2 ln 2 (2) -iln 4 (2) +i]Vln 4 (2) -12a 4 + ^iVa4 



(65) 



with £ mM = ln(/i 2 /M 2 ). 
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